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We study in a Brill-Hartle type of approximation the back-reaction of a superposition of 

linear gravitational waves in an Einstein-de Sitter background up to the second order in 

the small wave amplitudes hjfc. The wave amplitudes are assumed to form a homogeneous 

QJ ^ and isotropic stochastic process. No restriction for the wavelengths is assumed. The 

^ , effective stress-energy tensor T/;^ is calculated in terms of the correlation functions of 

^\ ■ the process. We discuss in particular a situation where T' is the dominant excitation of 

^^ ' the background metric. Apart from the Tolman radiation universe, a solution with the 

^\ ' scale factor of the de Sitter universe exists with p = —p as effective equation of state. 

T-H . 

^ ' While the description of cosmic gravitational radiation in linearized relativity 

04 I is fairly well known, a study of its nonlinear aspects is much harder. Numerical 

J^^ ' relativitya is one way, but there are also some analytical and semianalytical ap- 

proaches. Already in 1964 Brill and HartleQ proposed a scheme which takes the 
back-reaction of linear gravitational waves into account. As it stands, the Brill- 
[~^ I Hartle method is considered as a high-frequency approximation for gravitational 

On ■ radiation. Thus its application to early inflationary stages of the Universe is ques- 

^ \ tionable, since low-frequency modes may be present, which turn into high-frequency 

O^' modes (with wavelengths smaller than the temporary horizon) only at later time. 

5^ I As shortly discussed in this note, a slight modification of the Brill-Hartle approach 

bX)' can remove this shortcoming. The detailed calculations will be published elsewhere. 

The wave perturbations are considered as random variables, forming a stochastic 
process, which shares the symmetries of the background metric. The method is sim- 
ilar to the Monin-YaglomQ approach to statistical fluid dynamics, and is also used 
in optical coherence theoryO. For simplicity, an Einstein-de Sitter model is cho- 
sen as background metric. Tensor perturbations are added in synchronous gauge, 
assuming goo = — a^, 9oi = Oj 9ik = a^Sik -\- hik, as well as the gauge conditions 
hii = 0, hik,k = 0. a{r]) is the scale factor. The Einstein tensor may be expanded 
in powers of hik- Retaining terms up to second order and performing a stochastic 
(ensemble) average on the field equations, they split into linear wave equations for 
the hik and the back-reaction equations. The back-reaction equations relate the 
Einstein tensor for the background metric to the effective stress-energy tensor of 
the waves, which is represented as stochastic average over bilinear expressions in 
hik and derivatives of hik- As well known, the effective stress-energy tensor is not 
gauge invariant in general . However, as shown by Abramo, Brandenberger and 
MukhanovB (see also a, eI), the gauges change the background geometry to second 
order, and these changes just compensate the change of the stress-energy tensor. 
Representing the hik as stochastic Fourier integrals J7ifc(k, ri)e^^^dk-\-conj.compl., 
the amplitudes jik satisfy an ordinary differential equation. The symmetry proper- 
ties of the problem allow a simple representation of expressions which are bilinear 



in hik, e.g. 

{hikhlm) = — rr— (3(5iZ(5fcm + iSimSkl ~ 2SikSlm) / dkk^ f , (1) 

in terms of a single spectral density f{k,ri). f satisfies a nonlinear differential 
equation (eo depends only on fc = |k|, and a prime denotes the differentiation with 
respect to 77) 

2ff"-r'+Af{k'--)-Aeo^0. (2) 

a 

For high-frequency waves /cry S> f , / = y^/fc follows as solution ("high-frequency 
approximation" ) . It is convenient to write p and p in the effective stress-energy ten- 
sor in terms of four frequency-independent, but in general time-dependent integrals 
("moments") over the spectral density f{k,'q): 

fo^ Jdkk^^-^, f,^ Jdkk^^, f2 = Jdkk^f, U^ Jdkk^f. (3) 

For a general scale factor one obtains (g = /q + /i/4) 

3P. = ^(7A-5, + 5^/^-5^/.). (5) 

In the high-frequency approximation iGa^Pg — Ga^pg — J dkk^^/eQ^ energy density 
and pressure are positive. If low- frequency modes are present, their contribution can 
be negative, and also the equation of state can deviate considerably from the high- 
frequency relation p = p/i. If pressure and density of gravitational waves cannot be 
neglected compared to other forms of matter, the back-reaction on the scale factor 
must be taken into account. Taking a pure gravitational radiation universe, one has 
to solve the field equations 

(6) 
(7) 
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Note the further equations 


9' = 


a 


f^' - 


2.g + 2-/2 - 2/4, 
a 



with pg , pg taken from 

(8) 

(9) 

which follow from differentiating /i,/2 and using the differential equation for /. 
The four equations (^)-(^) give the relation 

(2a72 - a/2) {aa" - 2a'^) = 0. (10) 
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Vanishing of the first factor leads to the Tolman radiation universe, vanishing of 
the second factor gives the scale factor of the de Sitter cosmos. It is easy to find the 
moments f2,fi,g from the corresponding differential equations. Self-consistency 
however requires, that the moments found in this way must be compatible with 
the expressions following directly from (y|), if the solution of the wave equation 
is inserted. Compatibility can be achieved indeed, it however requires singular 
infrared components for some spectral quantities. The general solution of (g) for 
the radiation (s = 0, a ^ 77) and de Sitter (s = 1, a ~ l/*?) cosmos is known to be 

f = 2npp* + {I + im)p^ + {I - im)p*^ , (11) 

where I, m, n are three functions of k, connected with eg by eo = 4fc^(n^ — P — m^), 
p{x) = {l + is/x)e''^ is a complex function oix = krj. Whereas n{k) is not restricted, 
the spectral functions l{k) and m(fc) should be understood as generalized functionsEl 
and have the form {b is a constant)" 

,.,^ bS"{k) 6ik) 

m{k) = 0, (13) 

where S(k) is the Dirac delta function and n-z = J dkk'^n{k). The spectral densities 
for the energy density and the pressure are then (again in the Tolman case) 



a^Gp{k,r,) - n{k){k'-—^) + 5{k){^-b), (14) 

ia^Gp{k,7^) = n{k){e-^) + 5{k)(^~b). (15) 



These quantities show infrared singularities, part of the effective energy density 
and pressure resides in infrared {k = 0) modes. It is so far not clear whether the 
singularities result from using only a second-order approximation. The integrated 
quantities are finite, however. 
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for the Tolman cosmos, the condition is diflferent for the de Sitter case. 
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